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Recently, three reactor neutrino experiments, Daya Bay, Double Chooz and RENO have 
directly measured the neutrino mixing angle On. In this paper, another important oscillation 
parameter, effective Am^i (= Ani^i) is measured using baseline dependence of the reactor neutrino 
disappearances. A global fit is applied to publicly available data and Arh^i = 2.95!t]^o7 ^ 10~ ^^ ' 
sin^ 2^13 = 0.099looig are obtained by setting both parameters free. This result is complementary 
to Amli to be measured by spectrum shape analysis. The measured Am|i is consistent with Am|2 
measured by j/^ disappearance in MINOS, T2K and atmospheric neutrino experiments within errors. 
The minimum x^ is small, which means the results from the three reactor neutrino experiments are 
consistent with each other. 



I. 



INTRODUCTION 



Neutrino oscillation is a phenomenon which is not 
accounted for by the Standard Model of elementary 
particles, which assumes neutrinos as massless. There 
are six parameters in standard three flavor neutrino 
oscillation [ij]: three mixing angles between flavor 
eigenstates and mass eigenstates (6*12, 6*13 and 6123), one 
CP violating imaginary phase (S), and two independent 



In reactor-013 experiments, usually the neutrino 
disappearance is analyzed by a two flavor neutrino 
oscillation formula; 



P{ve -> i^e) = 1 - sin^ 2^13 sin' 



Am^ 



31 



4^;, 



L, 



(1) 



Am2^ 



mj. , where 



squared mass differences: 
rui are neutrino masses (mi, TO2, rn^) of the three 
mass eigenstates {vi, 1^2, 1^3) which correspond to 
the largest component of (z^e, Vp,, Vt), respectively. 
9i2 and ^ra\2 have been measured by solar neutrino 
disappearance experiments {ve -^ t'e) and long baseline 
reactor neutrino disappearance experiments {D^ -^ 
Ve)- ^23 and Am32 have been measured by {v^ -^ 
Vp) disappearance experiments at accelerators and 
atmospheric experiments. All these measurements 
are summarized in [ij . Here, Am^ is a weighted 
average of Am^^ and Am|2, called effective Am^ as 
described in detail later in this section. Recently, 
finite 013 was finally measured by short baseline reactor 
neutrino experiments j^^land long baseline accelerator 
experiments {v^j, -^ Ve) [al^- 

Another effective mass squared difference AtTi^]^ can 
be measured by energy spectrum distortion and baseline 
dependence of the reactor-0i3 experiments. This paper is 
to measure Amgj^ by baseline dependence of the reactor 
ncutrino-013 experiments. 



where L is baseline which is ~ 1 km and E^ is neutrino 
energy, which is around a few MeV. Amgj^ is a weighted 
average of the two mass square differences, |Ato|]^| and 
^32! 



IATO32I of the standard parametrization. 



Am3i 



C12 Am3J 



S12 ATO32 , 



(2) 



with Cij and s^- representing cos 9ij and sinOij, respec- 
tively 3. In the analyses of reactor-6'13 experiments so 
far published, sin^ 26*13 is extracted assuming Am^i = 
Am32, which is measured by MINOS experiment [8|. 
Am32 can be expressed as. 



Arh^2 = (*i2 + 513^23 sin 2^12 cos 6) Am^i 

+ (Ci2 - 313^23 sin 2012 COS S) I ATO32 I 

where tij = tan 9ij [7| . Since there is a relation 



(3) 

(4) 

in the standard three neutrino flavor scheme, the 
difference between Arh^i and Ato§2 is expressed as 
follows, 

2(Am|i - Ami 



Am^i — A?7i32 + Attiji 



31 32^_^(^_^^^^^^^^^j^20i2COS(5) 

2 cos 26*12 1 Ato^;^ I 



Am^j^ + Am|2 
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where the overall sign depends on mass hierarchy. The 
difference is much smaller than the current precisions of 



measurements and can be treated practically equivalent. 
However, if Am|j^ and Am^j are separately measured and 
if they turn out to be significantly different, it means the 
standard three flavour neutrino scheme is wrong. Thus it 
is important to measure Am|]^ independently from Am§2 
to test the standard three flavour neutrino oscillation. 

The E dependence and L dependence analyses to 
extract Am^^ use independent information, namely 
energy distortion and normalization and thus are 
complementary. Some of the authors demonstrated 
Aml^ measurement using L dependence of deficit value 
of each reactor-0i3 experiment in 2012 0, [lO|. In this 
paper the analysis is significantly improved by applying 
a detailed global fit making use of the publicly available 
information of the three reactor neutrino experiments. 

In next section we re-analyze the published data of 
each experiment and compare with the results written in 
the papers in order to demonstrate our analysis produces 
identical result. Section-Ill discusses about possible 
correlations between the experiments. In section-IV, 
most recent Double Chooz, Daya Bay and RENO 
results 0, |j, UM ^^re combined and Am^i is extracted. 
Finally, a summary of this study is presented in section- 
V. 



II. 



REACTOR NEUTRINO DATA 



Details of each experiment and their data are presented 
in this section and they are re-analysed by the authors in 
order to demonstrate that the analysis methods used in 
this work are consistent with the publications from the 
experimental groups. The x^ used in this section will be 
used to form global x^ function in section-IV. 
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FIG. 1. Relative locations of detectors and reactors of Daya 
Bay Experiment. Scale is approximate. 



where Md are the measured neutrino candidate events 
of the d-th AD with background subtracted, Bd is 
the corresponding background, Td is the prediction 
from neutrino flux, Monte Carlo simulation (MC) and 
neutrino oscillation. ujf is the fraction of neutrino 
event contribution of the r-th reactor to the d-th 
AD determined by baselines and reactor fluxes. The 
uncorrelated reactor uncertainty is cr^. Cd is the 
uncorrelated detection uncertainty, and as is the 
background uncertainty, with the corresponding pull- 
terms {arjedTVd)- An absolute normalization factor a 
is determined from the fit to the data. 

The values of cu'} are not shown in Daya Bay 
publications and was estimated using 



Pr/L] 



2 
rd 



J2riPr/Lld) 



with Pr — 



'Er^r' 



(7) 



A. Daya Bay 

The Daya Bay (DB) reactor neutrino experiment 
consists of three experimental halls (EH), containing 
one or more antineutrino detectors (AD). The AD array 
sees 6 reactors clustered into 3 pairs: Daya Bay (DBl, 
DB2), Ling Ao (LI, L2) and Ling Ao-II (L3, L4) power 
stations. Fig.-[T] shows the relative locations of reactors 
and AD and table-U shows the distance between each 
combination of reactor and detector. All reactors are 
functionally identical pressurized water reactors with 
maximum thermal power of 2.9 GW jSj]. 

In DB publication, the x^ is defined as 



XDB(^i3,Am^i) = 



E' 



Md 



rid-Td[l + a + J2r ^i'^r + <^rf 



H^TT^^T. 



Md + Bd 



vj 



(6) 



where Wr is the thermal power of each reactor and Lrd 
is the baseline of r-th reactor to d-th detector. In 
this analysis, the value of Pr is considered 1/6 since 
all reactors have same nominal thermal power. The 
calculated w^ is shown in table-|lll All the others terms 
are shown in table- lHll 

By using equation (|6]) and the data from tables-lHl and 
- lllli we were able to reproduce Daya Bay's result, where 
Td was multiplied by the value of the deficit probability 
(Pdf), defined as: 



P;^f = 1 - sin^ 261 



Ji.sMcV Sin (1-27 



^lllJ^)n,{E)dE 



Jn^{E)dE 



(8) 

with Arrr^ being measured in eV'^, Ldr in meters and 
E in MeV. n^(E) is the expected energy spectrum of 
the observed neutrinos which is calculated by ni,{E) = 
S{E)aiBD{E). S{E) is the energy spectrum of the reactor 
neutrinos, which is a sum of the energy spectrum of 
neutrinos from the four fissile elements: 



S{E,) 



E 



|3^S^{E,), 



(9) 



i— 235]J 238U 23gpy 241py 



where Si{E^) is reactor neutrino spectrum per fission 
from fissile element i and /?i is a fraction of fission rate 
of fissile element i. For equilibrium light water reactors, 
Pi are similar and we use the values of Bugey paper |12| , 
namely ^ssy : 238u . 239py . 24ip^ ^ q 533 . q g^g 

: 0.328 : 0.056. In this study, Si{E) is approximated as 
an exp onential of a polynomial function which is defined 
in 



13| 



Si{E^) oc exp 



E 



a,4^"-i) 



(10) 



CTiBD is the cross section of the inverse process of 
neutron /3-decay (IBD), that can be precisely calculated 
from the neutron lifetime [IJ]. The energy dependence 
of the IBD cross section is. 



<Jibd{E,) (X (^.[MeV] - 1.29)^E^-2.59E, + 1A. 

(11) 
sin^ 26*13 is extracted by fixing Arn^i as the MINOS 
Am§2 = 2.32 X lO'^ eV^ [8]. The x^ distributions of 
the Daya Bay paper and our calculation arc compared 
in fig. -[2] The Daya Bay central value and uncertainty is 
sin^ 2013DB — 0.089 ± 0.011 while our analysis showed 
sin^ 2^13 = 0.090l;QQi[|!), in good agreement with the 
published value. We also verified how different values for 
the fission rates coefficients of equation (fTO|) and different 
assumptions for equation-©, affect the final result. It 
does not change the result significantly. 
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FIG. 2. X distribution with respect to sin 2613 by fixing 
Arn^ as Aifi^2 for Daya Bay data. The black curve is the x^ 
distribution shown in their paper ll| with central value and 
la uncertainty of 0.089±0.011, while the red curve shows the 
X^ distribution calculated in this analysis with central value 
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RENO. Scale is approximate. 



Together with the distances of each detector reactor 
pair, the contribution of each reactor flux to each detector 
for the period of their first analysis is available [l^ and 
are summarized in table- ITVl 

The RENO x^ is defined as: 



XRE(^i3,Am2^) = 



E 



[^^ob. 



U)Eli^ + fr)N,% 
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E Jr 
rr2' 



(12) 



B. RENO 



where Nl!;^^^ is the number of observed IBD candidates 
in each detector after background subtraction and 
N^^p is the number of expected neutrino events, 
including detection efficiency, neutrino oscillations and 
contribution from the r— th reactor to each detector 
determined from baseline distances and reactor fluxes. A 
global normalization n is taken free and determined from 
the flt to the data. The uncorrelated reactor uncertainty 
is ar, (7^ is the uncorrelated detection uncertainty, and 
CT^ is the background uncertainty, and the corresponding 
pull parameters are {fr, ^d, bd). The values of these 
variables are shown in the tableJVl 

The expected number of events for both detectors are 
not present in the RENO paper, but the ratio between 
data and expectation is shown. This ratio was used to 
calculate the expectation value (iV^ ). 

Using the data in the table-|Vl equation-® and the 
MINOS Ato§2> we obtained sin^ 26I13 = 0.111 ± 0.024 
which is in good agreement with their published value 
of sin^ 26'i3R,E = 0.113 ± 0.023. The x^ distributions are 
also very similar as shown in fig.-|4l 



The Reactor Experiment for Neutrino Oscillation 
(RENO) is located in South Korea and has two identical 
detectors, one near (ND) and one far (FD) from an array 
of six commercial nuclear reactors, as shown in fig.-|31 



C. Double Chooz 

The Double Chooz (DC) experiment uses the two 
Chooz B reactors with thermal power of 4.25 GWth each. 




FIG. 4. x^ distribution with respect to sin'^ 2613 by fixing 
Am as Am32 for RENO data. The black curve shows the x 
distribution shown in their paper [J| with central value and 
1 a uncertainty of 0.113±0.023, while the red curve shows the 
X^ distribution calculated in this analysis with central value 
and 1 a uncertainty of 0.111 ± 0.024. 



Currently, the experiment is using only the far detector, 
since its near detector is not complete yet. The Bugey-4 
measurement [12] is used as a reference of the absolute 
neutrino flux in the analysis, and the relative location of 
the far detector and reactors are shown in fig. -[51 where 
the distances from the detector to each reactor are 998.1 
and 1114.6 meters llal. 
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Daya Bay and RENO cases, assuming that the detector 
and background related uncertainties of [13] and are 
fully correlated. 
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where N°^^ is the number of the observed neutrino 
event candidates. The subscript "d" represents reactor- 
ON and OFF period. A^*"^? is the number of expected 
neutrino events, including detection efficiency and 
oscillation effects, and B is the total expected number 
of background events. The ar, ad, and ab are 
the reactor, detection and background uncertainties, 
respectively. The corresponding pull parameters are 
(a, e, b). Using the parameters shown in table- TVTl we 
obtained sin^ 26*13 — 0.131 ± 0.048 which is consistent 
with the result of the DC publication, sin^ 26*13 = 0.109± 
0.039, although the background evaluation methods are 
different using different data sets. We also did a rate 
only analysis of the Double Chooz data, which result 
agreed with the published one. Thus the reproduction 
of the Double Chooz data is also well performed with 
our analysis method. 



FIG. 5. Relative locations of detector and reactors of Double 
Chooz experiment. Scale is approximate. 



III. CORRELATION EVALUATION OF 
SYSTEMATIC UNCERTAINTIES 



The Double Chooz collaboration published a rate 
plus shape analysis result [2]. An effect of the 
shape analysis in this case is an evaluation of main 
backgrounds of ^Li and fast-neutron from the energy 
spectrum beyond the reactor neutrino energy range. 
Since information of detailed energy spectrum, which 
is necessary to reproduce the analysis, are not publicly 
available, we do not consider here the shape analysis 
but restrict only to the rate analysis. After the 
second publication on the sin 2^13 measurement, the 
Double Chooz group published a result of the direct 
measurement of backgrounds by making use of 7.53 days 
reactor-OFF period [ij]. We used these data in addition 
to the background evaluation inputs written in (2| 
to improve the background estimation instead of the 
energy spectrum analysis. The relative neutrino-flux 
uncertainty for reactor-OFF period is much larger than 
reactor-ON period. The dominant uncertainty comes 
from long-life isotopes whose abundance are not well 
known. It has negligible contribution in reactor-ON 
period [17|. Therefore, we regard the error correlation on 
neutrino flux between the reactor ON and OFF periods 
to be uncorrelated. We performed a similar x^ analysis as 



In reactor neutrino experiments, the expected number 
of observed events (A^cxp) is defined by: 



^-P - 4^i2 



A^pe^W), 



(14) 



where L is the reactor-detector baseline, iVp is the 
number of targets in the detector, e is the detector 
efficiency, Pth is the reactor thermal power, {Ef) is the 
mean energy released per fission, and (a/) is the cross- 
section per fission defined as: 



W) 



Y,|3^Js,{E)amBiE)dE. 



(15) 



For each experiment, L, A'p, e, and Pth terms are deter- 
mined independently. Therefore they can be assumed to 
be uncorrelated. On the other hand, (Ef) and (ct/) terms 
are taken from the same references and the uncertainties 
of these terms are correlated between the experiments. 
From the Bugey and Chooz experimental results, the to- 
tal uncertainty on spectrum prediction is 2.7%, where a 
2% correlation is expected between the experiments as 



treated in [l8|. Fully correlated signal prediction uncer- 
tainties between experiments, which come from neutrino 
flux and detection efficiency, can be cancelled by overall 
normalization factors used in the analyses of the Daya 
Bay and RENO. It allows us only to take into account 
remaining uncertainties between detectors or periods for 
each experiment. Daya Bay and RENO treat the remain- 
ing uncertainties as uncorrelated in their publications. 



IV. COMBINED ANALYSIS 

As explained before, the main method of this work is 
to combine all the data of the current neutrino reactor 
experiments in a single x^ function. Then we look for the 
minimum x^ value, calculate the Ax^ distribution, and 
determine the confidence level regions. The y^ function 
used for such analysis was chosen so as to use the data 
from tables-Uto - I VII as well as the correlation as described 
in section- Hnl The definition of our global x^ is, 



In fig. -[8] the baseline dependence of the disappearance 
probability of each detector is shown, where the 
probability is calculated using the parameters output 



Xg = Xdb + Xdc + Xre> 



(16) 



with the y^ of each experiment defined as in section- 
II. Therefore, this function has 32 pull terms: 18 for 
Daya Bay (6 reactors, 6 detectors and 6 backgrounds), 
10 for RENO (6 reactors, 2 detectors and 2 backgrounds) 
and 4 for Double Chooz (2 reactors, 1 detector and 
1 background). It also contains the two overall 
normalization factors, one for Daya Bay and the other 
for RENO data set. 

For all combinations of Am^ and sin^ 26*, the Xq is 
minimized with respect to the pull terms. Fig. -15] shows 
map of the absolute x^ and fig. -El shows Ax^ contour 
map near the Xmim obtained by such procedure. From 
the minimum point and its 1 a error region, 

Xmin = 5.14/6 Degrees of Freedom, 
Am^, = 2.95+?:^? (+°;^4,) X 10-3 eV^ 
sin2 20i3 = O.O99+n?^(+°:°l^) 

are obtained. The errors without parentheses are the 
full range of the two dimensional 1 a confidence level 
region, set when Ax^ = 2.30, for a determination of two 
free parameters fit. The uncertainties in the parentheses 
are the 1 a range when the other parameter is fixed to 
the best fit value. This Am§]^ is consistent with Amgj 
measured by accelerator experiments [8|, Il9| , confirming 
the standard three fiavor neutrino oscillation within the 
error. The sin^ 26'i3 obtained here is independent from 
Amij. The small Xmin/I-'oF means the data from the 
three reactor neutrino experiments are consistent with 
each other. 

All the pull terms output were within 1 a from the 
input value, and the normalization factors obtained from 
the fit to the data, were both less than 1%. 
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FIG. 6. Absolute y^ value calculated for each pair of Amli 
and sin^ 1Q\z, and by the minimization of the pull terms. For 
higher values of Am3i (bigger than IQ-^eV^) some valleys 
are present, although they are about more than ten times less 
sensitive than the minimum x^- 




FIG. 7. Ax^ distribution around the Xmin- From the inner 
to the outer part, the lines correspond to 1 a, 2 a and 3 a 
confidence level. The star shows the best fit point. There 
is no solution more significant than 3 a except for the Xmin 
valley. 



which give the best fit. The Double Chooz has a large 
effect on this Amgj^ determination because it locates at a 
baseline where the slope of the oscillation is large. In the 
near future, when the near detector of the Double Chooz 
experiment starts operation, the accuracy of this Am|]^ 
measurement is expected to improve much. 
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FIG. 8. Reactor v^ survival probabilities. The solid line is the oscillation pattern obtained in this analysis and dot-dashed line 
uses MINOS Am,32 and the sin^ 20i3 that returns the minimum x^- The data points are below the Am32 because they are 
calculated using the parameters returned by the best fit solution. Generally, a detector sees several reactors. The horizontal 
axis is a weighted baseline (L) and the horizontal bar in each data point shows the standard deviation of the distribution of 



the baselines, which is defined by gi, = \J iL)'^ — (L^), where (i") 
and Pk are the baseline and thermal power of the reactor k. 



lluP^Ll 



I "Ylih PkL^. . k is the reactor index and Lk 



Complenientary to this study, we demonstrated a 
similar, but simpler and robust measurement of the 
effective Ato|]^ from the baseline dependence of the 
disappearance probabilities of the three reactor-^is 
experiments [9|, [lOj. The result obtained on that work 
of Am|i = 2.99ti:58 ^ 10"^ ^V^ is compatible with the 
value obtained in this paper. In addition, a similar Ax^ 
distribution is presented in [20, fig--4]. 
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SUMMARY 



In this work, a global fit of the data from all the current 
reactor-013 experiments was performed to measure 
Am|]^. The combination of the data from Daya Bay, 
RENO and Double Chooz resuhed in Am|i = 2.95l?;^?x 
10"'^ eV^. This is consistent with Am^2 ^^'^ it confirms 
that the experiments are observing standard three flavor 
neutrino oscillations within the error. The mixing angle 
obtained this analysis is sin^ 26*13 = O.Ogglgoig. The 
small Xmin/DoF value indicates that the data from the 
three reactor experiments are consistent with each other. 
This analysis uses independent information from the 
energy spectrum distortion and it is possible to improve 
the accuracy of Arfi^i combining with results from energy 
spectrum analysis. It will be impormant to perform 
this kind of analysis to improve Aifi^^ accuracy and to 
check the consistency of the results from the reactor-0i3 
experiments. 
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I. Daya 


Bay: 


Baselines, 


in meters 


, between each 
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TABLE II. Daya Bay: 


Contribul 


;ion to each detector from 


detector and core 


y. 










reactor 


;-.-. 














DBl 


DB2 


LI 


L2 


L3 


L4 




DBl 


DB2 


LI 


L2 


L3 


L4 


ADl 


362 


372 


903 


817 


1354 


1265 


ADl 


0.4069 


0.3854 


0.0654 


0.0799 


0.0291 


0.0333 


AD2 


358 


368 


903 


817 


1354 


1266 


AD2 


0.4089 


0.3870 


0.0643 


0.0785 


0.0286 


0.0327 


AD3 


1332 


1358 


468 


490 


558 


499 


AD3 


0.0330 


0.0318 


0.2676 


0.2441 


0.1882 


0.2354 


AD4 


1920 


1894 


1533 


1534 


1551 


1525 


AD4 


0.1208 


0.1241 


0.1894 


0.1892 


0.1851 


0.1914 


ADS 


1918 


1892 


1535 


1535 


1555 


1528 


AD5 


0.1201 


0.1248 


0.1895 


0.1895 


0.1847 


0.1913 


AD6 


1925 


1900 


1539 


1539 


1556 


1530 


AD6 


0.1209 


0.1241 


0.1892 


0.1892 


0.1851 


0.1914 







TABLE III. Daya 


Bay: Fitting paramet 


ers [H. 








ADl 


AD2 


AD3 


AD4 


AD5 


AD6 


V candidate 


69121 


69714 


66473 


9788 


9669 


9452 


3do 


68613 


69595 


66402 


9922.9 


9940.2 


9837.7 


BKG total/day 


13.68 ±1.54 


13.55 ±1.54 


10.38 ±1.17 


3.56 ±0.24 


3.55 ±0.24 


3.44 ± 0.24 


Live Time (days) 


127.5470 


127.5470 


127.3763 


126.2646 


126.2646 


126.2646 


Efficiency 


0.8015 


0.7986 


0.8364 


0.9555 


0.9552 


0.9547 


Ma 


67723.59 


68334.17 


65363.96 


9358.7 


9240.98 


9037.24 


(JB 


157.43 


156.86 


124.65 


28.95 


28.94 


28.93 


0-d 


0.002 


0.002 


0.002 


0.002 


0.002 


0.002 


Or 


0.008 


0.008 


0.008 


0.008 


0.008 


0.008 



TABLE IV. 
tions flal. 



RENO: Basehnes and neutrino flux contribu- 



Rl 



R2 



R3 R4 R5 R6 



FD baseline [m] 1556.5 1456.2 1395.9 1381.3 1413.8 1490.1 

FD contribution 0.1373 0.1574 0.1809 0.1856 0.1780 0.1608 

ND baseline [m] 667.9 451.8 304.8 336.1 513.9 739.1 

ND contribution 0.0678 0.1493 0.3419 0.2701 0.1150 0.0558 



TABLE V. RENO: Fitting parameters [i]. 



ND 



FD 



IBD candidates 


154088 


17102 


^ *cxp 


151723.54 


17565.72 


BKG (total/day) 


21.75 ±5.93 


4.24 ±0.75 


Live Time (days) 


192.42 


222.06 


Efficiency 


0.647 


0.745 


A^obs 


149902.86 


16160.46 


Od 


0.002 


0.002 


Or 


0.009 


0.009 


Ob 


1141.05 


166.54 



TABLE VI. Double Chooz: Fitting parameters [l[i3l- 



Reactor-On 



Reactor-Off 



IBD candidates 


8249 


8 


IBD prediction 


8439.6 


1.42 


BKG (total/day) 


2.18 ±0.58 


2.00 ±0.58 


Live Time (days) 


227.93 


6.84 


Od 


0.010 


0.010 


Or 


0.017 


0.40 
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